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Abstract 

We complete and improve the fits to experimental tttt scattering am- 
plitudes, both at low and high energies, that we performed in the previous 
papers of this series. We then verify that the corresponding amplitudes sat- 
isfy analyticity requirements, in the form of partial wave analyticity at low 
energies, forward dispersion relations (FDR) at all energies, and Roy equations 
helaviKK threshold; the first by construction, the last two, inside experimental 
errors. Then we repeat the fits including as constraints FDR and Roy equations. 
The ensuing central values of the various scattering amplitudes verify very ac- 
curately FDR and, especially, Roy equations, and change very little from what 
we found by just fitting data, with the exception of the D2 wave phase shift, for 
which one parameter moves by 1.5 a. These improved parametrizations there- 
fore provide a reliable representation of pion-pion amplitudes with which one 
can test various physical relations. We also present a list of low energy param- 
eters and other observables. In particular, we find 4°' = 0.223 ± 0.009 



-0.0444 ± 0.0045 and <5^°'(m^) - 5'^^\m\) = 50.9 ± 1.2°. 
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1. Introduction 

In two recent papers,[^'^l that we will consistently denote by PY05, KPY06, we have presented a 
set of fits to the data on tttt scattering phase shifts and inelasticities, and we also checked how well 
forward dispersion relations (FDR, henceforth) are satisfied by different tttt scattering phase shift analyses 
(including our own), finding in fact that some among the more widely used sets of phase shifts, failed 
to pass this test by several standard deviations. We then performed (in PY05) a consistent energy- 
dependent phase shift analysis of tttt scattering amplitudes in which we constrained the fits by requiring 
verification of FDR. 

In the present paper we complete and improve on the results in PY05 and KPY06 in various 
ways. First of all (Sect. 2) we incorporate the fit to data in ref. 3 (to be called GMPY07 henceforth) for 
the SO wave which is very precise thanks to use of recent experimental results^'*! and a more appropriate 
parametrization, which provides an accurate determination of 5q^^ (s) below KK threshold. The precision 
of this determination is such that we have to refine also the analysis of some other waves. Thus, 
we improve the S2 wave as given in PY05 by requiring smoother junction between the low energy 
(s^/^ < 932 MeV) and intermediate energy (932 MeV < s^/^ < 1420 MeV) regions. We also (slightly) 
improve the fit to the inelasticity of the DO wave, smoothing the onset of inelasticity above the KK 
threshold, and the parametrization of the P and F waes removing their ghosts. Finally, we improve the 
error analysis of the D2 wave. 

We also make (Sect. 3) some improvements in the high energy (s^/^ > 1.42 GeV) input, espe- 
cially in the momentum transfer dependence of the amplitudes, necessary to evaluate Roy equations. 
Indeed, for Roy equations, we need Regge formulas away from the forward direction, so we extend our 
analysis of the Regge amplitudes there. For the imaginary part of a scattering amplitude, hn F{s,t), 
one only expects Regge theory to be valid for s ^ A^, \t\ ^ s (with A ~ 0.35 GeV the QCD parameter); 
this will limit the validity of the Roy equations to energies ^ 1 GeV (as a matter of fact, we here 
only evaluate them below KK threshold). The corresponding Regge parameters were obtained, in the 
forward direction, by relating the tttt cross section to the ttN and A'^A'^ cross sections, using factorization, 
and fitting all three processes. Away from the forward direction, we fix other parameters using also sum 
rules, in Sects. 5 and 6. 

All of this (Regge parameters and partial wave amplitudes) constitute a set of pion-pion ampli- 
tudes, obtained with unconstrained fits to data, that we denote by "UFD Set" , which we have collected 
in Appendix A. 

Next, in Sect. 4, we test FDR and Roy equations for this UFD Set. We get good verification in 
both cases, to a level of average agreement corresponding to 1.00 standard deviations, for FDR. and to 
0.97 standard deviations for the Roy equations. The fact that FDR and Roy equations are practically 
satisfied within errors, makes it reasonable to improve the fits to data including as a constraint the 
verification of forward dispersion relations and Roy equations: this should provide a set of parametriza- 
tions to partial waves fully compatible with analyticity and, hence, with more reliable central values. 
This we do in Sect. 5; here we also constrain the fits by requiring verification of two sum rules which 
relate higher (s-'^/^ > 1420 MeV) and lower energies, and which permit us to refine the values for Regge 
parameters away from the forward direction. This procedure provides a set of tttt amplitudes, that we 
call constrained fit to data, or CFD Set, that not only fit data, but also satisfy sum rules, FDR and Roy 
equations, well within the rather small errors we now have; in particular, the degree of verification of the 
FDR below 932 MeV and, especially, of the Roy equations, is spectacular. It turns out that all param- 
eters for all waves remain inside the error bars we obtained by just fitting the data (UFD Set), except 
for the D2 wave where, in particular, one parameter changes by 1.5 cj. The result of these constrained 
fits is collected in Appendix B for the partial wave amplitudes and Regge parameters. As it happened 
in PY05, we still find a marked hump in the SO wave both for the unconstrained and constrained fits, 
between 400 and 900 MeV of energy, which structure is thus shown to be compatible with both FDR 
and Roy equations. In Sect. 6 we present the values of the low energy parameters and other observables 
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that follow from both fits to data (UFD Set) and, especially, constrained fits (CFD Set). The values of 
some parameters are improved using sum rule determinations. In particular, for the scattering lengths 
for the SO, S2 waves we find the very accurate determination 

4°' = 0.223 ± 0.009 M-\ = -0.0444 ± 0.0045 M'^; (1.1) 

here, and in what follows, represents the mass of the charged pion, that we take = 139.57 MeV . 
The article is finished in Sect. 7 where we give the conclusions, and a short discussion of our results; 
particularly, comparing them with other independent evaluations (theoretical and experimental) of the 
TTTT scattering amplitudes. A few words on isospin breaking corrections (which affect very little our 
results) for the SO, P waves are also said there. 

2. Fits to data 

In the present Section we briefly summarize the methods and results in PY05, KPY06, together with a 
few improvements, obtained fitting experimental data for tttt partial waves, up to the energy of 1.42 GeV. 
Above this energy we assume the scattering amplitudes to be given by Regge formulas, that we discuss 
in Sect. 3 below. The reason why we choose to use Regge theory above 1.42 GeV is that, whereas 
there are uncertainties in the Regge expressions at the lower part of the energy range where we use 
them, say between 1.42 GeV and ~ 1.8 GeV, these are substantially smaller than the uncertainties in 
TTTT phase shift analyses in the same energy region. Broadly speaking we distinguish two energy regions 
in the fits to experimental partial wave data: the low energy region, s^^'^ < sY'^, s|^^ ~ 1 GeV, and 
the intermediate energy region, < 5^/2 < GeV. is the energy above which one cannot 

1 /2 

consider that inelastic processes are negligible. The precise value s^' where we separate "low energy" 
from "intermediate energy" depends on each wave. 

In the low energy region we write model-independent parametrizations that take into account 
unitarity and analyticity; this last, by making a conformal mapping, 

/ s — \/ Si — s 

s wis) = —= , 

\/s + V-Si - s 

and expanding the function^ coi ?)f^ i^s) in powers of this variable w. The details may be found in 
Appendix C, where we also explain how using this variable does not imply any model, while accelerating 
the convergence: we, generally speaking, only need between two and four terms in this expansion. We 
then fit experimental tttt phase shift data^'*'^! and, for the SO wave at intermediate energy, also tttt — > KK 
data.I^l In four cases (the waves DO, D2, F, GO) we include in the fit the values of the scattering length 
and, for the DO wave, also the effective range parameter. These are obtained from experiment via the 
Froissart-Gribov representation; see PY05, KPY06 for details. In the case of the P wave, we do not 
fit the data on tttt scattering, but the vector form factor of the pion, which gives much more precise 
results. In the intermediate energy region we make phenomenological fits to experimental phase shifts 
and elasticity parameters, basically polynomial fits. There is, however, an exception to this: the SO wave. 
Here, most of the inelasticity is due to the KK channel, so we can make a two-channel calculation; see 
KPY06. The precision of the fits, particularly for the SO, P and DO waves, is such that we require a 
(slight) improvement of the analyses of PY05, KPY06 in some cases. 

^ Up to a kinematical factor, and up to poles, that we separate explicitly; see the specific formulas for each wave. 
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2.1. The SO wave 

The first case where we improve on the analysis of PY05, KPY06 is the SO wave, at low energy. 
Although this is discussed in detail in GMPY07, we say a few words here for completeness. To write the 
corresponding expansion, it is necessary to separate off the pole of the effective range function that lies 
on the real axis, viz., the pole due to the so-called Adler zero of the partial wave amplitude. This lies 
near the beginning of the left hand cut, at s = |2q, — -^tt- We fit using the expression 

We here fix zq — M^^; in Sect. 5 we will allow zq to vary. We find, taking the best fit^ in ref. 3, 

So = 4.3 ±0.3, Bi = -26.9 ± 0.6, S2 = -14.1 ± 1.4. (2.1b) 

The ensuing numerical results for the phase shift are very similar to what was obtained in PY05, 
but have a sounder theoretical basis and, above all, are much more accurate. The curves corresponding 
to this, and that in PY05, are shown in Fig. 2.1, together with some experimental data. Because at 
932 MeV we match the low energy to the K-matrix fit at intermediate energy, and the low energy fit 
has changed (in particular with much smaller errors) from what we had in KPY06, we have also to 
slightly modify the parameters of the K-matrix fit. The resulting values for the parameters are given 
in Appendix A. They do not change much from what we had in KPY06, but they are now determined 
more accurately. 

The corrections due to isospin breaking have been computed for this wave, at very low energy; 
we discuss the (slight) modification they imply in Sect. 7, and the corresponding values of the parameters 
are given also in Appendix A. 

2.2. The S2 wave 

The second case where we improve on the analysis in PY05, KPY06 is for the S2 wave. The reason is 
that the precision of our calculations has improved so much that one is sensitive to the derivative of the 
phase shift as one crosses into the region where inelasticity is not negligible. In the case of the S2 wave, 

the lowest inelastic process is tttt irnp, which is a (quasi-)three body process. Therefore, we expect 

(2) 

6q (s) to be continuous, and with a continuous derivative, until the first two-body channel [pp) opens. 
We take this into account now by altering the PY05 fits to this wave, as follows. In the low energy 
region, s^/-^ < 932 MeV, we maintain the fit in PY05. We write 

"'*i"W = ^ J^{^» + '''^T^}' = '•"'i G-Vi (2.2a) 

< 932 MeV . We fix Z2 — M^; later, in Sect. 5, we will allow Z2 to vary. For the errors, and since 
the low and high energy pieces are very strongly correlated, we take this into account and get somewhat 
improved errors. Altogether we find a /d.o.f. — 11.2/(21 — 2) and the parameters 

So = -80.4 ±2.8, Bi = -73.6 ± 10.5. (2.2b) 

^ In ref. 3 we gave several different parametrizations for the SO wave. These different parameters do not cor- 
respond to different physical scenarios, but just to the use of two or three terms in the conformal expansion, 
to what data sets are fitted, or to wheter or not we explicitly factorize the zeros in the amplitudes. Such 
parametrizations correspond to the same physical scenario, as indeed their resulting phase shifts overlap within 
errors. We here choose, among those parametrizations, the one with better analytic properties, and which fits 
best the largest sample of consistent and reliable data points in the elastic region. 
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Figure 2.1. The fits to SO from GMPY07 (solid line and gray band) and in 
PY05 (dotted lines). In the indent, a blow-up of the low energy region. The 
fitted data from refs. 4, 6 are also shown. 

For the high energy region we neglect the inelasticity below 1.45 GeV for the fit to the phase. 
We then fit high energy data (s^/^ > 0.95 GeV), requiring agreement of the central value and derivative 

1 /2 

with the low energy determination given in (2.2) at the energy Sj^ — 932 MeV. We write 



5^/2 > MeV; 



hi 



Bho —Bo + Biwi{sm), Bhi — Bi 



Wh{s) - Wh{sM)] + Bh2 [Wh{s) - Wh{sM)f^ 

dwi{s) 



(2.3a) 



dwh{s) 



and 



^,(,) = ^7^gZI ; s]'^ = 1050 MeV, 



Wh[S) 



\/s + ^Jsi - s ' 

\/s - \/Sh - S _ 1/2 



1450 MeV. 



Bh2 is a free parameter. We get a reasonable /d.o.f. = 13.8/(13 — 1) and 

Bh2 = 112 ± 38. 



(2.3b) 
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Figure 2.2. The fits to the S2 phase shift here (solid line and gray band) 
and in PY05 (dashed lines). Some experimental data from ref. 6 are also 
shown. 
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Both the present fit and that in PY05 may be found in Fig. 2.2. 

The present fit has slightly smaller errors than the old fit in PY05 in the low energy region, and 
larger ones above 1 GeV. 

The inelasticity we still describe by the empirical fit in PY05: 

1 - e(l - si/sf^, e = 0.18 ± 0.12; s > si = (1.05 GeV)^; 
1, s<si = (1.05 GeV)2. 



ryf (s) 



(2.4) 



2.3. The inelasticity for the DO wave 

We now consider the elasticity parameter for the DO wave. In PY05, KPY06 we made a phenomenological 
fit to r?f': 



1 - e 



s < 4m^, 
k2{s) 



■; s>4m|.; e = 0.262 ± 0.030, k2 = ^/s]T 



(2.5) 



This provides a fit to the elasticity parameter on the average. 
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Figure 2.3. Fit to rl^^ (continuous line and dark area that covers the 
uncertainty), from Eq. (2.6) here. The dotted lines follow from KPY06 [Eq. 
(2.5) here]. The elasticity on the /2(1270), from the Particle Data Tables, 
is also shown as a large white dot; other data are from ref. 6. 



The problem with (2.5) is that it rises too brusquely at KK threshold,^ proportional to /c2, 
which causes distortions for s^/^ near 2mx\ a behaviour /cf near this threshold would be indicated. For 
this reason, we here also try the following parametrization: 



s < 4m^, 



< 



1 



1 



5/2 



4m|- 



-5/2 



his) 



(2.6) 



s > 4m|.; 
=0.284 ± 0.030; 



m 



K- 



2.54 ±0.31; M 



f2 



1275 A MeV . 



This will underestimate the inelasticity near 2mK, since, in fact, most of the inelasticity of the DO wave is 
not due to KK but to the four pion states; but will have a smooth threshold behaviour. Moreover, (2.6) 
is more flexible than (2.5) and it will, by construction, have the correct inelasticity around the /2(1270) 
resonance, which is the most important region. If evaluating the dispersion relations with (2.5), (2.6) we 



^ We thank H. Leutwyler for this remark. This type of problem is important for the DO wave, because KK 
threshold is not far from the very strong /o(1270) resonance; but it also exists for P and D2 waves. We have 
checked that, in the case of these waves, the influence of the incorrect threshold behaviour is much smaller 
than our errors, hence negligible. 
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find that they are almost equivalent. We will use (2.6) in our calculations in the present paper because, 
as stated, it reproduces best the inelasticity around the /2(1270) resonance. 

The elasticity parameter here and that from PY05 are shown in Fig. 2.3. 



2.4. The D2 wave over the whole range 



Finally, we discuss the D2 wave. The data for the corresponding phase shift are very poor (and nonexis- 
tent for the inelasticity). For this reason, we made, in PY05, a phenomenological fit in the whole energy 
range from threshold to 1.42 GeV. We wrote 

cot 6i'\s) -'^{Bo + BMs) + BMsf} 4(^2^'^-. (2.7a) 
with A a free parameter fixing the zero of the phase shift and 



Wis) — 



= 1450 MeV, 



and we got a mediocre fit, /d.o.f. = 71/(25 — 4); the values of the parameters were 

Bo = (2.4±0.3) X 10^ = (7.8±0.8)xl0^ B2 = (23.7±3.8) x 10^ Z\ = 196± 20 MeV . (2.7b) 

When we fitted requiring fulfillment of dispersion relations, the parameters moved to (PY05, Ap- 
pendix A) 

Bo = (2.9±0.2) X 10^ Bi = (7.3±0.8) X 10^, B2 = (25.4±3.6) x 10^ Zi = 212± 19 MeV . (2.7c) 

The corresponding phase shift moves by a little more than one standard deviation. This shows that the 
errors in (2.7b) were underestimated. We improve this by keeping the central values of the parameters 
we found in the fit to data, but enlarging the errors by adding quadratically the difference between the 
central values in (2.7a) and (2.7b), which may be considered as a "systematic" error in the fit.'* Thus 
we find the parameters 

So = (2.4±0.5)xl0^ Bi = (7.8±1.0)xl0^ B2 = (23.7±4.2) x 10^, Zi = 196±25 MeV . (2.7d) 

It should be noted that the errors in ^2^^ are still small, compared to the experimental errors; (2.7d) 
gives errors for below the level of 0.8° . 

For the elasticity parameter we write the same formula as in KPY06: above 1.05 GeV, 

r^f'^ (s) = 1 - e(l - s/sf, s^/^ = 1.05 GeV, e = 0.2 ± 0.2. (2.7e) 



* It should be remarked that the present method, as indeed any method for estimating "systematic" errors, 
is arbitrary to a large extent. In PY05 we took into account the poor quality of the fit (which reflects the 
incompatibihty of the various sets of experimental data) by scaling the purely statistical errors by the square 
root of the x^/d.o.f., ^71/(25 - 4) ~ 1.8. This produced (2.7b). We could have scaled instead by the x^/d.o.f. 
itself; this would have given errors like in (2.7d) for the parameter, Bo, and a bit larger than those in (2.7d) 
for the other parameters. Likewise, we could have taken as central values the averages of the central values in 
(2.7b) and (2.7c). We consider our method to be reasonable and, as we find in the present paper, the increase 
in the errors quite justified. 
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2.5. The GO wave 



With respect to the GO wave, we gave in PY05 a parametrization for its imaginary part based on 
dominance by the resonance /j, 



k{Ml) j {s- Mlf + Mlr^[k{s)/k{Ml)]^^ ' (2.8) 
s^/^ > 1 GeV; BR = 0.17 ± 0.02, M/^ = 2025 ± 8 MeV, T = 194 ± 13 MeV . 

Unfortunately, this much underestimates the value at low energy. In fact, from the Froissart-Gribov 
representation we can evaluate rather accurately the scattering length, finding af^ = (8.0 ± 0.4) x 
10~^ M~^. At low energy one can write 

Im/f (s)2.[fc(s)]l«[af]^ (2.9) 

which disagrees with what one would find from (2.8a) by many orders of magnitude. A simple formula 
that interpolates between low and high energy is 



18 



I„/<».W= Jf) BR 



2 2rd-s/.A/f )- 



k{Ml) J {s- MlY + Mlr^[k{s)/k{Ml)r ' (2.10) 

BR = 0.17 ±0.02, M/^ = 2025 ± 8 MeV, r = 194±13MeV; c = 9.23 ±0.46. 

This gives the correct values near the /4 resonance and at low energy, the only regions where we have 
experimental information. Note that this interpolation is rather arbitrary, but there is no point in 
trying to improve it as there are no data to which one can fit more realistic formulas. Any reasonable 
interpolation would give the same order of magnitude estimate for the contribution of this wave to 
dispersion relations and sum rules. 

We have, in our calculations, not taken into account the contributions of GO or G2 waves other 
than to check that they are considerably smaller than the experimental errors: this is the only interest 
of the parametrization. 

2.6. Low energy SO, P and F waves: ghost removal 



When cutting the low energy expansions 

cot 6{s) = K{s){Bo + Biw + • • •} 

(with K an appropriate kinematical factor) at a finite order, a ghost, i.e., a spurious pole in the partial 
wave amplitude /(s) appears in the vicinity of the point s = for the SO, P and F waves. As remarked 
in GMPY07 (see especially Appendix A there), such ghost poles are rather harmless, their effect being 
at the percent level: removing the ghosts is little more than an aesthetical requirement. Nevertheless, 
we will here improve our formulas by writing them in such a manner that the ghosts disappear. This 
was already done for the SO wave in GMPY07, and the SO parametrization used in the present paper 
takes this into account; now we remove also the ghosts for P and F waves. As a matter of fact, this 
improves the consistency of our results, slightly but systematically. 

In PY05 and KPY06 we used the following formula for the P wave, 

cot<5i(s) = |^(M2-s)|i?o + Bi^^^^|; 4/^ = 1.05 GeV. (2.11a) 

The best result, from ref. 8 (see also PY05), is 
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So =1.069 ±0.011, Si = 0.13 ± 0.05, Mp = 773.6 ± 0.9 MeV. (2.11b) 
Instead of this we now write a parametrization where the ghost is absent, 

coti,(.) = l^(M,^-.){J|| + B„ + B,^^^}; 4'=^l«GeV, (2.12a) 

and find 

Bo = 1.055 ±0.011, Bi = 0.15 ±0.05; (2.12b) 

the difference with what follows from (2.11) is less than 0.7% at tttt threshold, decreasing to 0.05% at 
KK threshold. 

For the F wave, one can remove the ghost without changing the parameters we found in PY05 
(within the significant digits), so we have 

-«»)-^--{^--« + --^^^}' ^r = 1.45GeV; ^^^^^ 
So = (1.09 ± 0.03) X 10^ Si = (1.41 ± 0.04) X 10^ A = 0.051 x 10^ 

2.7. Other waves 

The results of the remaining fits for this UFD Set may be found in PY05 and KPY06, with details of 
the fitting procedure and the far from trivial matter of the selection of data. The results of all the fits 
are collected in Appendix A, for ease of reference. 

2.8. Matching 

Before turning to the calculations of forward dispersion relations, and Roy equations, a few words have 
to be said about the matching points between our low energy and intermediate energy regimes, and 
between intermediate energy and high energy regime, that we discuss later. Such matchings are, of 
course, artificial: for example, and as discussed about the DO wave, the inelasticity is small, but not 
zero, below KK threshold. A fully satisfactory matching is however not possible; it would require 
a multichannel evaluation, and hence introducing a number of parameters impossible to fix with the 
existing experimental data. In KPY06, and here, we have requested matching of the central values of 
the phase shifts, dl^\s), at the matching points themselves, usually (but not always) the KK threshold. 
However, and because the errors in the low energy and intermediate energy ranges are independent, this 
produces jumps (when varying the parameters inside their error bars) which increase the errors of the 
dispersive integrals artificially. This is, unfortunately, an unavoidable feature of our analysis: the results 
deteriorate somewhat when one is very near the matching points. 

In fact, the situation is even less clear near the KK threshold. In our analysis we neglect isospin 
breaking effects, and therefore we have taken it at an average between the K^K^ and K'^K'^ thresholds, 
si/2 = MeV. Since the K+K^ and K°K° thresholds differ by some 8 MeV, the threshold itself is 
thus not well defined to this extent: 992 ± 4 MeV. All in all, the net result is that our dispersive (or 
even direct calculations) of the tttt amplitudes suffer from uncontrollable errors in a, fortunately narrow, 
band of less than or around 6 MeV (to be on the safe side) around the matching points. We have avoided 
these matching regions when calculating the fulfillment of dispersive relations. 

With respect to the matching between intermediate and high energy regions, the situation is 
different. It is clear that, near s^/^ — 1420 MeV, which is the corresponding matching point, the Regge 
expression can only agree with the real amplitudes in the mean (as can be seen in the cases where we 
have precise data, as for ttN scattering). We expect that, since these Regge amplitudes only appear 
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in integrals, the fluctuations will be averaged out to lie within the errors. However, we here have a 
problem similar to that of the low energy matching with intermediate energy: this lack of correlation of 
the errors causes artificially enhanced error bars near the matching points. Therefore, our calculation 
should be used excluding a (narrow) band below the matching point. Here we have refrained to compare 
calculations above 1400 MeV, which is sufficient to render most of the fluctuations smaller than the 
experimental errors. 

3. Regge formulas 

Regge formulas have been obtained for tttt scattering, in the forward direction, by fitting experimental 
data for the various tttt total cross sections. This provides expressions that are not very precise. One 
improves this by use of factorization. It is then possible to include information on total cross sections for 
ttN and NN scattering!^! (A'' A'' includes antinucleon-nucleon scattering), which furnishes us with precise 
results for the contributions to tttt scattering of the three Regge poles^ P, P' and p. Here, we will use 
these Regge expressions above 1.42 GeV. 

These results are sufficient to calculate forward dispersion relations. For Roy equations, however, 
we require also the imaginary parts of the scattering amplitudes, Im F(s,t), for rather large values of 
\t\; in our calculation, up to t = —0.43 GeV^. In fact, these values are so large that one does not expect 
Regge theory to hold in the extreme range. What we do to circumvent this problem is to enlarge the 
errors in the t dependence of the parametrizations so that they cover, in the whole t range, all fits to 
experimental data. This will provide, at large t, a phenomenological representation of the corresponding 
scattering amplitude. For isospin zero exchange, we take the expressions given in PY05, in the forward 
direction, and enlarge the errors away from the forward direction. For p exchange we take, at t = 0, the 
parameters described in KPY06 and, for t 7^ 0, an uncertainty that covers the extreme fits in ref. 10 for 
vrA'' scattering, and assume that tttt scattering will vary in a similar manner. 



We write 



ImF(^*=i)(s,i) ^^^P, 




/3p =1.22 ±0.14, ap(0) =0.46 ±0.02; = 0.90 GeV'^; a^' = -0.3 GeV 
^{t) ^l + dpt + ept"^; 6 = 2.4±0.2 GeV-^. 



-4 



(3.1a) 



s is a scale parameter, that we consistently take s = 1 GeV^. We set 

dp = 2.4 ±0.5, ep = 0±2.5 GeV-''. 
For the Pomeron and P', one can write, also for s^/^ > 1.42 GeV, 

lmFii^='Hs,t) ^^^P{s,t)+P'{s,t), 



(3.1b) 



t fixed 




(3.2a) 



/3p = 2.54 ± 0.04, = 0.20 ±0.10 GeV"^ Wp{t) = I + cpt; 
pp, = 0.83 ± 0.05, ap, (0) = 0.54 ± 0.02, a'p, = 0.90 GeV"^; 



^p,{t) ^l + cpd; 6 = 2.4 ±0.2 GeV 



^ The P' is in fact a combination of two Regge trajectories, associated with the /2(1270) and /2(1525) resonances. 
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We may fix 

cp = 0.0 ± 1.0 GeV-^; cp/ = -0.4 ± 0.4 GeV"^ . (3.2b) 

If we do so, we cover the fits of Rarita et al.[^°l and of Proggatt and Petersen. I^^l 

Note that we do not give errors for the slopes a'p and a'p, because the variation of iZ'p/ covers 
possible variations of the Regge slopes: ^ varies a lot at large t. In fact, we have a range of variation 

-0.56 < ^(-0.4 GeV^) < 0.64, 

and something similar for Wp/. Fortunately, however, the evaluations for the Roy equations below 1 GeV 
depend very little on the scattering amplitudes for large s and large \t\. 
Finally, for exchange of isospin two we write 

ImF(^*=2)(s,t) ~ /32e''*(s/s)"''W+"''(°)-\ /32 = 0.2 ± 0.2; s > (1.42 GeV)^. (3.3) 

s—>oo 

These fits are expected to represent experimental data for energies between 1.42 GeV and 
^ 20 GeV and for 4M^ > t > —0.4 GeV^, with less reliability at the more negative values of t. At 
values of s^/^ larger than 20 GeV, one would have to use more complicated formulas, taking into account 
in particular the logarithmic growth of the total cross sections. For our purposes the formulas given 
above are sufficiently accurate, since the influence of the energy region much above 20 GeV for forward 
dispersion relations or Roy equations is negligible. 

The values of the Regge parameters can be improved by requiring verification of dispersion 
relations, and of two sum rules that relate directly the Regge behaviour, for nonzero t, to low energy 
amplitudes [see below Eqs. (5.1) and (5.2)]. The resulting numbers are collected in Appendix B. 



4. Dispersion relations 



4.1. Forward dispersion relations 



In this Section we will evaluate forward dispersion relations for the three independent tttt scattering 
amplitudes. For these calculations we will take the parameters for all partial waves from the fits to data 
described in the previous Sections (and collected in Appendix A). 

Although the form of the dispersion relations has been given before, we repeat them here. For 
7r°7r° scattering we write 

ReFoo(.)-Foo(4Mj- P-P-^/^^^ds 4^2) + 4m2) " ^^.la) 

The result of the calculation is shown in Fig. 4.1a, where the continuous curve is the real part evaluated 
from the parametrizations, and the dashed line is the result of the dispersive integral, i.e., the right hand 
side of (4.1a). 

The dispersion relation for n^ir'^ scattering reads, with fb+(s) the forward tt^tt'^ amplitude. 

In Fig. 4.1b we show the fulfillment of (4.1b). 

Finally, the dispersion relation for the It — ^ scattering amplitude does not require subtractions, 
and reads 

Re (. 0) - VP r ds' ^"^^'''"^(^^'Q) (4 ic) 

The result is shown graphically in Fig. 4.1c. 
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1- ReF„,(s)-ReF„,(4M;) 



s _1 - 




— direct 
dispersive 



Unconstrained Fits to Data (UFD) 
400' 600 800 ' 1000 




800 
s"^ (MeV) 



1200 1400 



Figure 4.1a. The tt^tt" dispersion 
relation. Continuous line: left hand 
side of {4.1a), evaluated directly with 
the parametrizations (the gray band 
covers the error). Dashed line: the 
result of the dispersive integral. 



ReF„^(s)-ReFJ4M^) 



-1-| ; 

— direct 
" --■ dispersive 

-2 -Unconstrained Fits to Data\ 
(UFD) 



400 600 800 1000 1200 1400 
s"^ (MeV) 



Figure 4.1.B The 7r°7r+ disper- 
sion relation with the new P and D2 
waves. Continuous line: left hand 
side of (4.1b), evaluated directly with 
the parametrizations. The gray band 
covers the errors. 

Dashed line: the result of the disper- 
sive integral. 




Unconstrained Fits to Data (UFD) 



400 



600 



8(30 
s"^ (MeV) 



1000 1200 1400 



Figure 4.1c. The dispersion 

relation for the /( = 1 amplitude. 
Continuous line: real part and error 
(shaded area) evaluated directly with 
the parametrizations. Dashed line: 
the result of the dispersive integral. 
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To measure quantitatively the fulfillment of the dispersion relations we evaluate the average 
(squared) distance between the real parts, calculated with our parametrizations, and the same real parts 
but now calculated with the aid of the dispersive integrals, a quantity that we denote by cP. This 
quantity is defined as follows. First, we rewrite the dispersion relations as the discrepancies Ai (shown 
graphically in Fig. 4.2), 



Ao+{s) ^ ReFo+(s) - Fo+i^M^) - ^ ^P.P. / 



s{s-4M^)^^ r (2,s'-4M2)ImFo+(s') 



Um^ s'(s' - s){s' - 4M2)(s' + s - 4M2) ' 

(4.2b) 

and 

These quantities would vanish, Ai = 0, if the dispersion relations were exactly satisfied. Because our 
fits have errors, we can only require vanishing within the uncertainties that said errors induce in the Ai, 
that we call SAi . Therefore, we define the quantities (average discrepancies) 

2 



* number of points ^ \SAi{sn) J 

This we do for all three relations (4.2). The values of the s„ are taken at energy intervals of 25 MeV. For 
the dispersion relation for It — 1, we also include the value at threshold, known at times as the (first) 
Olsson sum rule. For the other two dispersion relations, since the A vanish identically at threshold, we 
include a point below threshold, at s = 2M^. This is useful, among other things, to fix the location of 
the Adler zeros for the SO, S2 waves. 

If we had a fit to FDRs (which we do not) instead of an evaluation, (P would be the average 
chi-squared of the fit: in our case, (P is simply a measure of how well the forward dispersion relations are 
satisfied by the data fits, which are independent for each wave, and independent of dispersion relations. 
When calculating this (P , we use, in the present Section, the parameters for phase shifts and inelasticities 
discussed in the previous Sections, and collected in Appendix A. 

The average discrepancies in the various cases are given in Eq. (4.4) below: 

< 932 MeV 5^2 < U20 MeV 

7r°7r° FDR d? = 0.12 cP = 0.29 (4 4) 

77077+ FDR = 0.84 d^ = 0.86 

It = 1 FDR (P = 0.66 (P = 1.87. 

Below 932 MeV, the average (P shows a remarkably good fulfillment of FDRs. Still, the situation over 
the whole range is such that the (P for the It — i FDR is well above unity, indicating that there is room 
for improvement. 



AiiSn) 
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-0.5 - — Unconstrained fits to data (UFDy| 



400 600 800 1000 1200 1400 

1/2 

s (MeV) 




400 600 800 1000 1200 1400 
s'° (MeV) 



Figure 4.2. Fulfillment of dispersion relations, with the central parameters in (4.1a). The error 
bands are also shown. 
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4.2. Roy equations 

Roy equations are fully equivalent to nonforward dispersion relations, plus some t — s crossing symmetry, 
projected on the various partial waves. They can be written as 

r-OC , 

Re//^)(s) = CP4°) + C;(')42)+VP.P. / ds'Ki,i,,jj,is',s)lmfP{s'). (4.5a) 

Ci^\ C'l^^^ are known constants, and the kernels Ki^ii.jji are also known. 
We also define the quantities 

A\'\s) = RefPis) - - Cra^^^ -T^P-P- / ds' Ki,i,,jj,{s' , s)lm f^/ \s'y, (4.5b) 

they would vanish if Roy equations were exactly fulfilled. 

Roy equations are only valid up to s = 64 1 GeV^ because, for larger values, the integrand 

receives new contributions from the double spectral functions, not contained in (4.5). Moreover, when 
the value of t over which one integrates to project the partial waves is \t\ ^ ^qcd' for s ~ 2 GeV^, 
the Regge expressions are not valid. This is a further limitation of the validity of the Roy equations to 
energies below ^ 1 GeV. 

We here calculate up to KK threshold, and only test the waves SO, S2, P.*^ We define the 
equivalent of the average discrepancies (P we used for the forward dispersion relations, 

^ 1 ( ^Pi^n) 

''^ number of points ^ \6A['\sn) 

and find the results shown in Figs. 4.3, where we plot what our parametrizations give for Re f^^' (denoted 
by "in") and what follows from the integrals in the right hand side of (4.5a), denoted by "out". 
Numerically, we have the results 

4o = 0.54, 42 = 1-63, d|. = 0.74 : (4.7) 

a reasonable fulfillment, but with it is clear that this can be improved, particularly for the S2 wave for 
which the discrepancy is larger than unity. 

® A preliminary review of these results was presented at the 4th Int'l Conf. on Quarks and Nuclear Physics, 
Madrid, June 2006.1^31 
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Figure 4.3a. Fulfillment 

of the Roy equation for the SO 
wave. 

Continuous line: the result of 
the dispersive integral. 
Dashed line: real part, with 
the dark band the error band. 



Figure 4.3b. Fulfillment 

of the Roy equation for the S2 
wave. 

Continuous line: the result of 
the dispersive integral. 
Dashed line: real part, with 
the dark band the error band. 



Figure 4.3c. Fulfillment 

of the Roy equation for the P 
wave. 

Continuous line: the result of 
the dispersive integral. 
Dashed line: real part, with 
the dark band the error band. 
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5. Improvement of the parametrizations of the partial waves: 

fits to data constrained requiring fulfillment of dispersion relations 

5.1. Two sum rules 

Apart from forward dispersion relations and Roy equations, we will also require fulfillment, within errors, 
of two sum rules that relate high energy (Regge) parameters for t ^ to low energy P and D waves. 
The first sum rule is (PY05) 



1 = 



(5.1) 

The contributions of the S waves cancel in (5.1), so only the P, D, F and G waves contribute. At high 
energy the integrals are dominated by rho exchange. 

The second sum rule we consider is that given in Eqs. (B.6), (B.7) of ref 14. It reads, 



J AMI y 



; 4ImF^(°)(g,0)-10ImF^(^)(a,0) , ^ ImF^W(g, 0) - Imf W(s, 0) | 
/ s-{s-AMlY 6(3.-4mJ ,2(,_ 4^2)3 1 = 0. 

(5.2) 

Here F'^^\s,t) = dF^^\s,t)/d cos 6. At high energy, the integral is dominated by isospin zero Regge 
trajectories. We also define a discrepancy for these sum rules: 

These two sum rules are reasonably satisfied, if using the partial wave parameters obtained from data 
(Sect. 2), and the Regge parameters determined from factorization and fits to data (Sect. 3). 

5.2. Minimization procedure 

Because forward dispersion relations and Roy equations are satisfied almost within the fluctuations 
induced by the experimental errors, it makes sense to repeat the fits to experiment requiring verification 
within errors of forward dispersion relations and Roy equations, to which we add the sum rules (5.1), 
(5.2) to control Regge parameters away from the forward direction. We do this by minimizing the 
quantity defined, with self-evident notation, by 

Here, p^^^ are the parameters that we have found in the unconstrained fits to experimental data, and 
5pi are their errors. Thus the sum over pi runs over BnS, zeros {zq, Z2, ^), inelasticity parameters e„, r, 
etc, and over the K-matrix parameters for the SO wave. The presence of the sum ^^[{p^ — pT^)I^PiY^ 
of course, ensures fit to experimental data. 

The quantity W in (5.4) is a weight, which can be estimated in two different manners. First, it 
will serve to give each of the dispersion relations (FDR or Roy) a weight appropriate to the information 
that they carry. For example, for the FDR for scattering, the quantity ReFoo(s) — Foo(4M^) 

can be fixed giving the slope at s = 4M^, the value of ReFoo(s) — -foo(4-M^) and its derivative at 
each of the points s where its changes direction; and the same at the end point, s — 1.42 GeV^ 
(see Fig. 4.1.A): altogether, 13 values. Any reasonably smooth function that fits these 13 values is 
sure to follow ReFoo(s) - Foo{AMl) in ah the range: putting extra weight would be imposing redundant 
constraints. For other dispersion relations the number is a bit smaller; in general, a number 6 < < 13 
is obtained. An alternate method to find W is to increase it so that all dispersion relations are satisfied 
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within errors, that is to say, all the corresponding cP are less than or equal to unity. This occurs for 
W ~ 9. In our calculation we have taken W — 9, although we have verified that results practically 
indistinguishable are obtained ioT 7 < W < 12. 



5.3. Results for the constrained fits (CFD) 



The results of the fits to data, constrained by requiring fulfillment of the FDR and Roy equations plus 
the sum rules (5.1) and (5.2), that we call Set CFD, are summarized in Appendix B. This CFD Set is 
obtained by minimizing as defined in Eq. (5.4). In general, the parameters hardly change with respect 
to what we had from fits to data, Appendix A; but there are a few waves for which there are noteworthy 
alterations. First of all, we have the SO and S2 waves. Because now we are requiring verification of 
the FDR below threshold, we can leave the location of the Adler zeros zq, Z2 free. This produces some 
changes in the parameters B^, since there is a strong correlation between them and the location of the 
Adler zeros. However, the phase shifts themselves are practically identical to what we had in fits to 
data. Appendix A. Then we have the D2 wave, the only one that changes appreciably. It moves by a 
bit more than one standard deviation. The S2 wave also moves appreciably after constraining its fit; we 
discuss the two later on (Sects. 5.4 and 5.5). All other waves change so little that the difference between 
Sets UDF and CFD is almost inappreciable. 

For the present CFD Set, FDRs and Roy equations are, of course, better satisfied than before. 
For the FDR we find 



(5.5) 



sV2 < MeV s^/^ < ^420 MeV 

7r"7rO FDR (P = 0.13 (P = 0.31 

7r07r+ FDR = 0.83 = 0.85 

It ^ 1 FDR (P ^ 0.13 (P ^ 0.70, 
and, for the Roy equations, 

d|o = 0.23, (Ps2 = 0.25, d^p = 2x 10'^. (5.6) 

Furthermore, for the sum rules (5.1) and (5.2) we have — 0.02 and — 0.55. The verification of 
FDR and Roy equations are shown graphically in Figs. 5.1, 5.2 and 5.3. 

The overall average P is now substantially smaller than unity; particularly for the Roy equations. 
However, for the FDR for tt'^tt^ scattering is near unity. This is probably due to the D2 wave, which is 
likely still a bit away from its correct location, and to the P wave at energies above 1380 MeV, where 
our parametrization fails to take account of the p(1450) resonance. This indicates that one is at the 
limit of accuracy for experimentally-based parametrizations of scattering amplitudes. 

Besides sum rules and dispersion relations, another independent test of our amplitudes is the 
Adler sum rule that relates the pion decay constant to pion-pion scattering amplitudes with one pion 
off its mass sell. This has been recently evaluated^^^^ with our scattering UFD amplitudes, and a very 
satisfactory fulfillment of the sum rule is found; the discrepancy A,^ that measures the accuracy with 
which the sum rule is fulfilled (and which should vanish if it was satisfied exactly) is found to be 
= 0.021 ± 0.053. 
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Figure 5.1a. The 7r°7r° dispersion 
relation with the CFD amplitudes. 
Continuous line: real part, evaluated 
directly with the parametrizations (the 
gray band covers the error). 
Dashed line: the result of the disper- 
sive integral. 
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Figure 5.1b The -k^t^^ dispersion 
relation with the CFD amplitudes. 
Continuous line: real part, evaluated 
directly with the parametrizations. 
The gray band covers the error. 
Dashed line: the result of the disper- 
sive integral. 
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Figure 5.1c. The dispersion re- 
lation for It = 1 scattering with the 
CFD amplitudes. 

Continuous line: real part and error 
(shaded area) evaluated directly with 
the parametrizations. 
Dashed line: the result of the disper- 
sive integral. 
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Figure 5.2. Fulfillment of dispersion relations, with the central parameters in (4.1a). The error 
bands are also shown. 
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Figure 5.3a. Fulfillment 

of the Roy equation for the SO 
wave. 

Continuous line: the result of 
the dispersive integral. 
Dashed line: real part, with 
the dark band the error band. 



Figure 5.3b. Fulfillment 

of the Roy equation for the S2 
wave. 

Continuous line: the result of 
the dispersive integral. 
Dashed line: real part, with 
the dark band the error band. 



Figure 5.3c. Fulfillment 

of the Roy equation for the P 
wave. 

Continuous line: the result of 
the dispersive integral. 
Dashed line: real part, with 
the dark band the error band. 
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5.4. Results: comparison of UFD and CFD Sets 



We here present the comparison of our Sets UFD and CFD, i.e., what we have by directly fitting data, 
and what is obtained constraining the fits by imposing also FDRs and Roy equations.^ Note that in the 
CFD Set we have here only altered the central values; we leave the errors that follow from fits to data, 
i.e., we assume errors as in Set UFD. Note also that, in the following formulas, the parameters are as 

defined in Appendices A and B. 

SO wave. We have, for SO below 932 MeV, 

UFD, Zo = CFD, zq free 

50 4.3 ± 0.3 4.41 ± 0.3 

51 -26.7 ±0.6 -26.25 ±0.6 (5-7a) 
B2 -14.1 ±1.4 -15.8 ±1.4 

Zo M„ 166.1 ±4.2 MeV. 

Above 932 MeV, 



S2 wave. We now find, 



For the inelasticity. 





UFD 


CFD 




0.843 ± 0.017 


0.843 ±0.017 


a2 


0.20 ±0.06 


0.20 ±0.06 


Pi 


1.02 ±0.02 


1.02 ±0.02 


^2 


1.33 ±0.013 


1.33 ±0.013 


711 


3.10 ±0.11 


3.10 ±0.11 


712 


1.82 ±0.05 


1.81 ±0.05 


722 


-7.00 ± 0.04 


-7.00 ± 0.04 


Ml 


888 ± 4 MeV 


888 ±4 


M2 


1327 ± 4 MeV 


1327 ± 4. 



UFD CFD 
0.17 ±0.12 0.18 ±0.12 



(5.7b) 



UFD, Zo = CFD, Z2 free 

50 -80.4 ±2.8 -80.2 ±2.8 

51 -73.6 ±10.5 -69.4 ±10.5 (5-8a) 
Bh2 109 ± 38 120 ± 38 

Z2 145.0 ± 3.6 MeV . 



(5.8b) 



P wave. In this case we have kept the value of the p resonance mass fixed when imposing dispersion 
relations; thus, for both Sets UFD and CFD, Mp — 773.6 ± 0.9 MeV. For the remaining parameters 
below 992 MeV we find. 



UFD CFD 

50 1.055 ±0.011 1.052 ±0.011 (5-9a) 

51 0.15 ±0.05 0.17 ±0.05. 

Above 992 MeV, 



For the SO wave, we give here only the values obtained neglecting isospin breaking; the values of the parameters 
obtained taking isospin breaking into account may be found in Appendix B. 



- 22 - 



-THE PION-PION SCATTERING AMPLITUDE. Ill: IMPROVING THE ANALYSIS - 





UFD 


CFD 


Ai 


1.57 ±0.18 


1.50 ±0.18 


A2 


-1.96 ±0.49 


-1.97 ±0.49 


ei 


0.10 ±0.06 


0.09 ± 0.06 




0.11 ±0.11 


0.12 ±0.11 


of the resonance fixed at Mf^ — 1275. 




UFD 


CFD 


Bo 


12.47 ±0.12 


12.48 ± 0.12 


Bi 


10.12 ±0.16 


10.12 ± 0.16. 




UFD 


CFD 


Bhi 


43.7 ± 1.8 


43.5 ± 1.8 


€1 


0.284 ±0.030 


0.283 ± 0.030 


r 


2.54 ±0.31 


2.53 ±0.31. 



(5.9b) 



DO wave. We here keep the ma 
CFD. We have, below 992 MeV, 



(5.10a) 

Bi 10.12 ±0.16 10.12 ±0.16. 

Above 992 MeV, 

UFD CFD 

(5.10b) 

D2 wave. This is the only wave that changes substantially; see Fig. 5.4. We find now, 

(5.11a) 



For the inelasticity parameter. 



F wave. This wave is unchanged within our precision: 





UFD 


CFD 


Bo 


(2.4 ±0.5) X 10^ 


(3.1 ±0.5) X 103 


Bi 


(7.8 ± 1.0) X 10^ 


(7.9 ±1.0) X 10^ 


B2 


(23.7 ±4.2) X 10^ 


(24.7 ±4.2) X 10^ 


A 


196 ± 25 MeV 


205 ± 25 MeV . 




UFD 


CFD 




e 0.2 ±0.2 


0.15 ±0.2. 



(5.11b) 



UFD CFD 

Bq (1.09 ± 0.03) X 10^ (1.09 ± 0.03) X 10^ (^-l^) 
Bi (1.41 ± 0.04) X 10^ (1.41 ± 0.04) x 10^. 

Regge parameters. We only give the values of the Regge parameters that we have allowed to vary. The 
parameters correspond to the formulas in Sect. 3. 
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FlGURE 5.4. The / = 2, D wave phase shift. Dashed line: fit to data with Eq. (2.7d). Contin- 
uous line and shaded area: after improving with dispersion relations. Also shown are data points 
from ref. 6. 



UFD CFD 

CP (0.0± 1.0) GeV"^ (0.53 ± 1.0) GeV"^ 

cp, -0.4±0.4 GeV"2 -0.38 ± 0.4 GeV"^ 

I3p> 0.83 ±0.05 0.83 ±0.05 (5.13) 

ap'(O) 0.54 ±0.02 0.54 ±0.02 

Pp 1.22 ±0.14 1.30 ±0.14 



a 



(0) 0.46 ±0.02 0.46 ±0.02 



(32 0.20 ± 0.2 0.22 ±0.2. 

The only parameters that change appreciably (but both by only ~ 0.5 a) are cp, which was to be 
expected, and (3 p. 

5.5. Comments 

From Eqs. (5.7) through (5.13) we see that the changes in most waves induced by the constraints 
given by FDR and Roy equations are very small; in many cases, minute or even nonexistent within the 
accuracy of our formulas. There are a few exceptions. First of all, the BnS for the SO and S2 waves 
change because now the location of the Adler zeros is left free (although the phase shifts themselves 
move very little). Secondly, the Bq and Bi parameters of the P wave vary by 0.3 and 0.4 a respectively, 
and, at high energy, the parameter Ai changes by 0.3 a. And thirdly, the only wave that suffers changes 
by more than one sigma is the D2 wave, as was to be expected; the parameter, Bq, moves by 1.5 a. All 
the other parameters of the SO, S2 and P waves, as well as all the parameters of DO, F waves, change 
below the limit of relevance. 

Altogether, the stability of the fit against imposing FDR and Roy equations is remarkable, 
showing its robustness. The stability is not obtained at the cost of large errors; quite the contrary. 
Except for the S2 wave at intermediate energies and for the D2 wave, where the errors are larger, the 
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errors in the other waves are as smaU (for the P, F waves) or much smaller (by a factor ^ 3) than what 
we had in previous fits, in PY05, KPY06, even after improving with FDR. 

Finally, a few words may be said on the Regge parameters. The parameters for exchange of 
isospin zero are almost unchanged when requiring fit to FDR, Roy equations and sum rules. This was 
to be expected; they are very well determined from ttN and NN amplitudes using factorization. For 
exchange of isospin unity, only the parameter Pp changes appreciably, and this by ~ 0.5 a. This shows 
the high degree of compatibility between our amplitudes above and below 1420 MeV. 



6. Low energy parameters and other observables 



6.1. General 

We present in the following Table 1 the low energy parameters (scattering lengths and effective range 
parameters, in units of the charged pion mass) that follow from our calculations. Besides scattering 
lengths and effective range parameters, defined as 



where fj^^ is the partial wave of definite isospin / and angular momentum / (in the elastic region) and 
Sl^^ stands for its corresponding phase shift, we give the quantities ag""* — a'^^ and S^\'rn'^) — 5g^''(m^). 
These quantities are relevant for pionic atom decays and kaon decays. We also give the combination 
2a[,°^ - 5a[,^^ that appears in the Olsson sum rule. 



2a, -5a, - SM^ J ^^^^ ds —^—^^ , (6.1) 
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which corresponds to the It — ^ amplitude forward dispersion relation evaluated at threshold. This sum 

(2) 

rule is very useful in determining a precise value for Oq ' . 

Besides the Olsson sum rule, the results from direct fits may be improved with the help of the 
Eroissart-Gribov representation (whose explicit form may be found in PY05), and the following sum 
rules: the sum rule [Eq. (6.8)] in PY05, which we repeat here in Eq. (6.2) for ease of reference, 



bi 



2 

1 



ds i 



_(S-4M2)3 s3 



ImF(^*=o)(s,0) + i 



(S-4M2)3 s3 



ImF(^*=i)(s,0) 



(s - 4M2)3 s3 



ImF(^*=^)(s,0) 



(6.2) 

and two sum rules involving the effective range parameters for the SO, S2 waves. These are obtained 
evaluating the limit as s ^ 4M2 of the ratio 

F(s,0) -F(4Af2,0) 
s{s - 4M2) ' 

in the forward dispersion relations for tt^tt" and 7r°7r+ scattering. We find, 

= , (2s'-4M2)ImFoo(sO 
,4M2 * s'(s' + s-4M2)(s'-4Af2)(s'-s)' 



b^°^ + 2b^Q^ = 6M^ lim P.P 



3af ' + 6^,^^ = 4M^ lim P.P. j 

S'-^M^ J 4 



ds' 



(2s'-4M2)ImFo+(s') 



(6.3) 



s'(s' + s-4M2)(s'-4M2)(s'-s)- 



Note that the limit has to be taken for s larger than 4M2 ; for s < 4M2 , the derivative of Re fo (s) diverges 
like (Constant)/!^. The principal part of the integrals is essential; the r.h.s. in (6.3) is convergent at the 
lower limit of integration only because it is a principal part, and one thus has 

P.P. / dx = 0, for y>0. 

Jo {x-y)vx 

Taking the value ai — (38.2 ± 1.3) x 10~3 from the Eroissart-Gribov representation, with 
CED waves (cf. Table 1 below) this gives 



b^°^ = 0.289 ± 0.008 M'^, ^^"^ = -0.081 ± 0.0035 M; 



(2) 



(6.4) 



6.2. The scattering lengths Oq°'', Oq^"* 



The results reported under the headings UFD, CFD in Table 1 are what is found by fitting experimental 
data on partial wave amplitudes. However, for the scattering lengths Cq"' and a'^^ one can improve 

the results using the Olsson sum rule. One takes the value of the combination 2ao°' — Soq^^ from the 
integral in (6.1), a value that is practically independent of that obtained fitting data, and thus sets the 
constraints 

,(0) 



24°^ - 5a'^' = 0.667 ± 0.018 [Olsson sum rule] 



a^o' = 0.223 ± 0.010 [CDE] 

= - 0.0451 ± 0.0088 [CDE] 



(6.5) 



(in units of M^r) 
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Unc. fits (UFD) 


Constrained (CFD) 


Sum rules, with CFD 


Best values 


"o 


0.231 ± 0.009 


0.223 ±0.010 




0.223 ± 0.009 




-0.052 ±0.010 


-0.0451 ± 0.0088 




-0.0444 ± 0.0045 **=' 




0.282 ± 0.014 


0.268 ± 0.014 




0.267 ± 0.009*-'^' 


24°) - 54^) 


0.720 ± 0.055 


0.672 ± 0.048 


0.667 ±0.018("> 


0.668 ±0.017 




51.7 ± 1.2° 


50.9 ± 1.2° 




50.9 ± 1.2° 





0.288 ± 0.009 


0.291 ±0.009 


0.289 ± 0.008*'*' 


0.290 ± 0.006 





-0.085 ± 0.010 


-0.084 ±0.010 


-0.081 ± 0.0035*^"*' 


-0.081 ± 0.003 


01 (x 10^) 


37.3 ± 1.2 


38.0 ± 1.2 


38.2 ± 1.3'*' 


38.1 ± 0.9 


61 (x 10^) 


5.18 ±0.23 


5.09 ±0.25 


5.42 ±0.91*''' 
5.13 ±0.19*°' 


5.12 ±0.15 


4°) (xlO^) 


18.7 ±0.4 


18.7 ±0.4 


18.33 ± 0.36*^' 


18.33 ± 0.36 


4'' (xlO^) 


2.5 ± 1.1 


2.4 ±0.9 


2.46 ± 0.25*''' 


2.46 ±0.25 


6f (xlO^) 


-4.2 ±0.3 


-4.2 ±0.3 


-3.82 ± 0.25*''' 


-3.82 ±0.25 


6(^) (xlO^) 


-2.7 ± 1.0 


-2.5 ±0.8 


-3.64 ± 0.18*''' 


-3.59 ±0.18 


as (x 10^) 


5.2 ± 1.3 


5.2 ± 1.3 


6.05 ± 0.29*''' 


6.05 ± 0.29 


ba (x lO'^) 


-4.7 ±2.6 


-4.8 ±2.7 


-4.40 ± 0.36*''' 


-4.41 ± 0.36 



Units of M„- The numbers in the "Sum rules" column are from the Olsson sum rule*"', 
the Froissart-Gribov representation,*''' the sum rule in Eq. (6.2),*°' and the sum rules in 
Eq. (6.3).*'" 

*°' This number is obtained composing the CFD values for Oq^' with the best value 
for 2ao°' — 54^^ ; which best value is obtained from the Olsson sum rule. *-'^' This number 
takes into account the "best values" for ttp^' , Oq^' given above. 

For the best values for ai, 61, 60°' and bg^' we have averaged what comes from constrained 
fits (CFD), with what one finds from the sum rules (since they are practically independent). 
The best values for the other parameters are as follows: from the constrained fits (CFD), 
for 4°'- For DO, D2 and F waves, the best values are those coming from the Froissart- 
Gribov representation; because our fits to data impose these values, it would not make 
sense to average them. However, we have averaged the results for 62 and 63 since we did 
not impose their values when fitting the data. 
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-0.035 






-0.04 






-0.045 






-0.05 






-0.055 






0.205 0.21 0.215 0.22 0.225 0.23 0.235 0.24 



Figure 6.1. The ellipses in the — plane corresponding to 1-sigma (thick continuous line) 
and 2-sigma (broken line). 



We can then fit the % , cIq minimizing (6.5). The resulting errors are strongly correlated (as 
is well known); the corresponding ellipse is shown in Fig. 6.1. We can uncorrelate the errors by using 
two new variables, x, y with 

4"^ =0.223 + 0.129 a; + 0.335 j/ 

= - 0.0444 - 0.335 X + 0.129?/; (6-6) 
x=0± 0.0087, 2/ = 0±0.027. 
This gives the central values, and errors, reported under the heading "Best values" in Table 1, 

4°^ = 0.223 ± 0.009, a[,^' = -0.0444 ± 0.0045. (6.7) 

This represents a reasonable improvement on the errors we had before for . We consider (6.7) to be 
our best result for these scattering lengths. 

The S-waves scattering lengths can also be compared with other experimental information, not 
used in our fits, that give directly the combination a^^^ — ag^'. Indeed, the value found in Table 1 for 
— o'q^ agrees very well with the following independent experimental determinations: from pionic 
atoms, [^^1 that give 

a[,") - = 0.280 ± 0.013 (St.) ± 0.008 (Syst.) M^^ 
and from K^,^ decays that implyl^^l 

- = 0.268 ± 0.010 (St.) ± 0.013 (Syst.) M^^ 
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Figure 7.1a. SO wave phase shift (CFD Set). Some data from refs. 4, 6 are also shown. Notice 
the hump around 700 MeV. 

7. Conclusions 

In the previous Sections we have given a representation of the tttt scattering amplitudes obtained fitting 
experimental data below 1.42 GeV, supplemented by standard Regge formulas above this energy, what 
we have called the UFD Set. We have shown that this UFD Set satisfies very well forward dispersion 
relations and Roy equations, as well as crossing sum rules. Then, we have improved the central values 
of our fits requiring, besides fit to data, verification of FDR, Roy equations and sum rules, getting what 
we have called CFD Set. The central values in this CFD Set lie well inside those of the UFD Set, except 
for the D2 wave. FDR are now very well satisfied, while the verification of Roy equations is spectacular. 

Altogether, we have confirmed the findings of PY05, KPY06 and (for the SO wave) of GMPY07: 
but we have now errors much smaller than in PY05. At present, the low energy P wave is known at the 
limit of validity of our formalism (a limit given by isospin breaking effects,^ that our analysis neglects); 
while the SO and DO wave are near this same limit. The SO, P and DO phase shifts, with the values of 
the CFD Set, are shown in Figs. 7.1a, b, c. 

In connection with the SO wave, there is some interest attached to the location of the so-called 
"(J pole". This has been discussed in GMPY07; we there found, with the parameters of the UFD Set, 

= 496 ± 6 (St.) ± 11 (Sys.) MeV, r^/2 = 258 ± 8 (St.) ± 2 (Sys.) MeV . 
^ Some isospin breaking effects will be discussed below. 
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Figure 7.1b. P wave phase shift (CFD Set); errors below 1 GeV are as the thickness of the hne. 
Some data from refs. 4, 6 are also shown. Note, however, that 5i is obtained fitting the tttt scatter- 
ing data only above KK threshold; below it, it is obtained fitting the pion form factor (cf. ref. 8). 



For the CFD Set we now have 

= 473 ± 6 (St.) ± 11 (Sys.) MeV, r„/2 = 257 ± 5 (St.) ± 2 (Sys.) MeV . 

Here the first error is the statistical one, due to the errors in the parameters of the SO wave, and the 
second the error induced by the extrapolation, as estimated in GMPY07. We gave the numbers not 
taking into account isospin breaking; if we took it into account (see below for a discussion of this), the 
central values would change a little, to Ma- = 478 and r'o-/2 = 242 for the CFD Set. 

Nevertheless, we must remark that a really precise determination of the location of this resonance 
requires the use of the Roy equations. We will present this in a future paper. 

The S2 wave is not found with a precision similar to that of the SO or P waves. This is, 
of course, due to the lack of accuracy of the experimental data on states with isospin two (but, on 
the other hand, the corresponding scattering length, , is found with good accuracy). The lack of 
experimental accuracy is what prevents more precise values for the D2 wave, and something similar 
happens for the F wave. 

By comparison with an analysis similar to ours, our values for SO, D waves are a factor of three 
or four more precise than those in ref. 14. Our D2 wave, uncertain as it is, is still much more reliable 
than what is given in ref. 14: where a parametrization for this wave is given that is incompatible with 
their own findings at low energy (the scattering length), and with Regge theory and experimental data 
at high energy. Finally, at high energy (in the Regge region), our amplitudes fit well the existing data, 
something that the amplitudes in ref. 14 fail to do by a factor of two. 
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Figure 7.1c. do wave phase shift (CFD Set); errors as the thickness of the hne. Some data from 
refs. 4, 6 are also shown. 

We have also used these scattering ampHtudes to evaluate low energy parameters for P, DO, 
D2, and F waves, and other observables, clearly improving on previous work. These parameters may 
then be used to test chiral perturbation theory to one and two loops, or to find quantities relevant for 
CP violating kaon decays. In particular, for the S-wave scattering lengths we find the very accurate 
determinations 

4°^ = 0.223 ±0.009 a[,^' = -0.0444 ± 0.0045 M^^ (7.1) 

This can be compared with what is found in ref. 18, using chiral perturbation theory with Roy 
equations and experimental data input, 

= 0.220 ± 0.005 M-\ a[,^' = -0.0444 ± 0.0010 M'K (7.2) 

that is, a remarkable agreement at very low energies. However, at higher energies our SO wave phase 
shift deviates somewhat from that obtained in ref. 18. This deviation is due to the hump we find in 
the 400 to 900 MeV region, that makes our phase shifts larger than those of ref. 18. [It is to be noted, 
on the other hand, that this hump is also generated naturally within the framework of unitarized chiral 
perturbation theory, ref. 19, when fitting all existing scattering data (using a large systematic error to 
cover all phase shift sets)]. The S2 wave also devites slightly (but significantly) from that of ref. 18 above 
500 MeV. 

One can ask the question whether it would be possible to improve on our precision. The answer 
is, no in the sense that our amplitudes agree, within errors, with theoretical requirements and with data.^ 

^ One may think that imposing chiral perturbation theory could lead to decreasing the errors of the tttt scattering 
amplitudes. However, this matter is compHcated and is left for a future publication: at least the analysis can 
be now made from the well grounded set of pion-pion amplitudes given by our fits. 
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Moreover, for some waves (notably, SO, P and DO, shown in Figs. 7.1) our precision is at the hmit of 
the estimated corrections due to breaking of isospin invariance. For other waves, known less precisely, a 
sizable improvement would require substantially improved experimental data; certainly for the S2 wave 
at low energy, and for all the waves above ~ 1 GeV. This is particularly true for the inelasticities, 
very poorly determined from experiment, except for the SO wave where the lucky coincidence that it is 
mostly given by KK states and the existence of data on the process tttt KK saves the day. It may be 
assumed that, perhaps, imposing exact fulfillment of FDR and Roy equations could improve our errors. 
We have found this impossible: increasing the weight of FDR and Roy equations in the CFD, Eq. (5.4), 
results on amplitudes that indeed satisfy better (but not that much better) FDR and Roy equations, 
but that deviate from experiment by intolerable amounts. This is not surprising; there are important 
features that fail to be resolved by existing experimental analyses. For example, we have the matter of 
the inelasticity of the DO wave near 1 GeV, or the lack of information on the /9(1450) resonance, clearly 
seen in e'^e" annihilations but for which the Particle Data Tables refrain from giving the 27r branching 
ratio -and which is absent from our analysis. This could give a sizable contribution above 1380 MeV: the 
absence of the /9(1450) resonance in our P wave is likely responsible for the structure found in Zio+i 
above ~ 1 GeV; cf. Figs. 4.2, 5.2. Unfortunately, taking this /9(1450) resonance into account correctly 
requires a multi-channel calculation, which, even if possible, lies outside the scope of the present article. 

As for the SO wave, the contribution of An states, dominant above ~ 1350 MeV is not properly 
taken into account. Likewise, our Regge formulas only give the imaginary parts of the scattering ampli- 
tudes in the mean, in the energy region 1420 MeV < < isoOMeV; a region where neither phase 
shift analyses nor Regge fits give precise results: our scattering amplitudes are not well determined in 
the region from 1350 to 1800 MeV. Until much better experimental data, in particular on phase shifts 
and (above all) inelasticities above KK thresholds are available, the Set CFD of pion-pion amplitudes 
will remain the best that one can do, from experiment. 

Next, we say a few words with respect to isospin breaking corrections, neglected in our analysis 
up till now. In most cases, our errors are sufficiently large to cover the estimated values of such effects; 
but there are two exceptions. For P wave, our errors include isospin breaking corrections. This is 
possible because we can evaluate this wave either from the pion form factor in +e~ — > 7r+ + tt^ or in 
r ^ 1/ -h 77+ -h 7r° : see ref. 8 for the details. And a special case is the SO wave. Here a recent calculation 
has been made^^^l in which account is taken, for Ke4 decays, of the fact that in the real world isospin is 
broken. According to ref. 20, this is done by subtracting, from the experimental phase shift, as given in 
refs. 4, the correction 

where i^V — 92.6 MeV is the pion decay constant, M^^o is the mass of the neutral pion, (3(s) — 
VI - 4M2/S, (3o{s) = ^l-U'P„/s and r = (m^ - m„)/(m, - md/2 - m„/2). 

If we take the results of this calculation at face value, we can repeat our fits taking it into 
account. For the UFD Set, only the low energy SO wave is affected, and it is so very little: by less 
than 1° above 400 MeV [which is not surprising, since, in fact, the correction (7.3) is actually somewhat 
smaller than the experimental errors over most of the range]. Only the scattering length moves beyond 
one standard deviation, to Oq^' = 0.210 ±0.010. For the CFD Set, the changes induced by incorporation 
of isospin breaking corrections in the waves other than the SO wave are almost negligible. The SO wave 
moves closer to the one given in the present paper, from which it is almost indistinguishable (Fig. 7.2); 

^°That the poor information above KK threshold is responsible for the limitations in the fulfillment of dispersion 
relations is seen if we realize that FDR are satisfied to an average dpDR = 0-49 below 932 MeV, and Roy 
equations are satisfied to an average cl^oy = 0-17 below KK threshold. 
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Figure 7.2. SO wave phase shift (CFD Set): continuous Une, neglecting isospin corrections; bro- 
ken hne, with isospin corrections taken into account. Notice that the hump around 700 MeV is ba- 
sically identical for both. Some data from refs. 4, 6 are also shown. 



and the scattering length becomes compatible with what we found neglecting isospin breaking: we now 
get 4°' = 0.213 ±0.011. Indeed, including isospin breaking corrections does not much affect our results. 
The corresponding values of the parameters for the SO wave, the only ones that change appreciably, are 
given in Appendices A and B; a comparison between what one finds with/without isospin corrections 
for the SO phase shift can be seen, with difficulty (because they are so close one to the other) in Fig. 7.2. 

We finish with a comment regarding the relative merits of Sets UFD and CFD. It is clear that the 
central values of the CFD Set should be considered as the preferred ones since they incorporate fulfillment 
of analyticity requirements: analyticity of the individual waves themselves (for the low energy region), as 
well as analyticity and crossing symmetry of the scattering amplitudes in the form of forward dispersion 
relations, sum rules, and Roy equations. However, the price to pay for the last two requirements is that 
all waves are now correlated: in this sense, the UFD Set is more robust. Since both solutions fit data, 
and are very similar, one can use one or the other almost interchangeably, except for S2 and, above all, 
D2 waves, for which the CFD Set is clearly superior. 
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Appendix A: fits to data up to 1.42 GeV (UFD Set) 



In this Appendix we collect the best values for the parametrizations of the various partial waves obtained 
by fitting experimental data, with the procedures defined in PY05 and KPY06; see also Sect. 2 here. 

A.l. The SO wave 

The region s^/^ < 932 GeV. We take sq — Am\, and take the Adler zero &t s — ^z^, with zq fixed at 
M^. We find 



cot4°'(s) = 



5I/2 Ml 



2k 



2^0 



Bq + Bi- 



Bo 



T 2' 



zo = M^; 



Bo =4.3 ±0.3, Bi = -26.7 ±0.6, B2 = -14.1 ± 1.4; 

(A.l. la) 

This corresponds to neglect of isospin breaking. If isospin breaking is taken into account as in 
ref. 20, we find instead 



,(0) 



Bo = 3.80 ±0.34, 
0.211 ± 0.010 M, 



-1 



Bi 



-27.1 ±0.8, 



B2 



0.278 ± 0.010 M, 



-3. 



(A.l.lb) 



The SO wave between 932 MeV and 1420 MeV. We here use the K-matrix fit of ref. 2." 
fci|fc2| det K ± kiKii 



tan (5^°^ (s) 



l + \k2\K22 

1 



2A;i[ifii +fc|if22 detK 



s < 4m|^, 

klKf:^ - klKl^ + klklideiKf 



(klKl^ + klKl^ + fc2A;2(detK)2 + 1)2 - Akjk^Kf, 



12 f 1 



s > AmJi 



and 



'{l + kik2 detK)2 + (fcij^ii -fc2X22)' 

(1 - kik2 det K)2 + {kiKu + k2K22Y' 



s > 4m%;. 



Here ki — ^ \/s^^TM^, k2 — \ \J s — 4m^ and the K-matrix elements are 



M2 



(A.1.3a) 



(A.1.3b) 



M{- s M2 - s fj, 

/i is a mass scale, that we take /x = 1 GeV. The powers of n have been arranged so that the a^, Pi, ^ij 
are dimensionless; they are also assumed to be constant. The pole at simulates the left hand cut of 
K, and the pole at A'/| is connected with the phase shift crossing 270° around 1.3 GeV; both poles are 
necessary to get a good fit. The values of the parameters are 



ai =0.843 ±0.017, 02 = 0.20 ± 0.06, /3i = 1.02 ± 0.02; /32 = 1.33 ± 0.013, 
711 = 3.10 ± 0.11, 712 = 1.82 ± 0.05, 722 = -7.00 ± 0.04; 
Ml = 0.888 ± 0.004 GeV, M2 = 1.327 ± 0.004 GeV; 



(A.1.3c) 



Ml lies near the beginning of the left hand cut for KK tttt scattering, located at 0.952 GeV. The 
parameters in (A.l. 3c) are strongly correlated. In fact, we have verified that there exists a wide set of 



^^A polynomial fit, giving results very similar to the K-matrix fit, can also be given. It may be found in 
Appendix B of KPY06. 
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minima, with very different values of the parameters. Nevertheless, the corresponding values of 6^^^ and 
77g°^ vary very little in all these minima, so that (A. 1.3c) can be considered a faithful representation of 
the SO wave for tttt scattering. The representations of the SO wave are matched exactly at 932 MeV, 
where one has 

4"^ ((932 MeV)2) = 104.9 ± 0.5° . 
A. 2. Parametrization of the S2 wave 

The region s^/^ < 992 GeV. For isospin 2, there is no low energy resonance, but this wave presents the 
feature that a zero is expected (and, indeed, confirmed by the fits). It is related to the so-called Adler 
zeros; to lowest order in chiral perturbation theory, one has the zero at s = 2z|, with Z2 = We 
note that, unlike the corresponding zero for the SO wave, 2z| is inside the region where the conformal 
expansion is expected to converge well. We here fix Z2 = and write 

cot^f (5) = ^ + , Z2 ^ M^; si'/' = 1.05 GeV. (A.2.1a) 
2k s - 2z^ [ ^/s + sjsi - s J 

Then we get 

So = - 80.4 ± 2.8, Si = -73.6 ± 10.5. (A.2.1b) 

The S2 wave between 932 MeV and 1420 MeV. We require junction with the low energy phase shift, 
and its derivative, at 932 MeV, neglect inelasticity below 1.45 GeV, and write 

cot (5^ '{s) = ^ _ |-Bfto + Bhi [wh{s) - Wh{sM)] + Bh2 [wh{s) - Wh{sM)] I , 

> MeV; (A.2.2a) 

dwi{s) 



Bho —Bo + Biwi{sm), Bhi — Bi 



dwh{s) 



wi{s) = ^ , ^{fjlZ ; SI = (1050 MeVf, 
y s + -s/ si — s 

WHis) - v^-^^fiH ; = (1450 MeV)2. 
Bh2 is a free parameter. We get 

Bh2 = 112 ± 38. (A.2.2b) 
The inelasticity is described by the empirical fit 

r^'^^\s)^l-e{l- si/sf/'^, e = 0.17±0.12 (s//^ = 1.05 GeV). (A.2.2c) 
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A.3. The P wave 

The region s^/^ < 2mK- We have 



,2_.J 2M^^;,„^«,^/W^I . 1/2 



cot.,(.) = ^(M; - .) 1^ + + ; . 1.05 GeV. (A.S.la) 

The best result is 

Bo =1-055 ±0.011, Bi = 0.15 ±0.05, Mp = 773.6 ± 0.9 MeV. (A.S.lb) 
The P wave between 2mK and 1420 MeV. We use a purely phenomenological parametrization: 

= Ao + Ai( 1) + A2(a/V4^- 1)', 

^ ^ (A.3.2a) 



77i(s) = 1 - €1^1 - 4m|./s - €2(1 - 4m|f/s); s > 4m%. 

The phase at the low energy edge, ((0.992 GeV)^) = 153.63 ± 0.55° , is obtained from the fit at low 
energy above; this fixes Aq. The rest of the parameters follow from the fit at intermediate energy. We 
have, 

Ao =2.681 ±0.010, Ai = 1.57 ±0.18, A2 = -1.96 ± 0.49; 

(A. 3.2b) 

ei =0.10 ±0.06, 62 = 0.11 ±0.11. ^ ^ 

A. 4. Parametrization of the DO wave 

The region s^^^ < 2mK- To take into account the analyticity structure, we fit with different expressions 
for energies below and above KK threshold, requiring however exact matching at s = Am\. Below KK 
threshold we take into account the existence of nonnegligible inelasticity above 1.05 GeV, which is near 
the pTTTT threshold, by choosing a conformal variable w appropriate to a plane cut for s > (1.05 GeV)'^. 
So we write 

1/2 

cot 4°) (s) = ^ {Ml -s)M^[Bo + B,w}, s< 4ml; 

(A.4.1a) 

w — —= , Sq' =1.05 GeV. 

V-s ± v-so - s 

The mass of the /2 we fix at Mf^ — 1275.4 MeV; no error is taken for this quantity, since it is negligibly 
small (1.2 MeV) when compared with the other errors. We find the values of the parameters 

So = 12.47 ± 0.12; Bi = 10.12 ± 0.16. (A.4.1b) 

The DO wave between 2mK and 1420 MeV. Above KK threshold we use the following formula for the 
phase shift: 

1/2 

cot4°'(s) = ^ {Ml - s) Ml{BhQ + Bh^w], s > 4m^; 

^ (A.4.2a) 

w = —p= , ; s, = 1.45 GeV . 

Vs± Vs^^ 

This neglects inelasticity below 1.45 GeV, which is approximately the pp threshold; inelasticity will be 
added by hand, see below. We require exact matching with the low energy expression, which yields the 
value of Bho- B^i follows from the fit at intermediate energy. We find 

Bm = 18.77 ± 0.16; Bhi = 43.7 ± 1.8. (A.4.2b) 
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For the inelasticity we write, as discussed in the main text, 
1, 



s < 4m|^, 



1 



4m|. 



5/2 



e =0.284 ± 0.030; r = 2.54 ±0.31. 



-5/2 



1 



k2{s) 
k2{Ml 



(A.4.2c) 



A. 5. Parametrization of the D2 wave 

For isospin equal 2, there are no resonances in the D wave. If we want a parametrization that applies 
down to threshold, we must incorporate the zero of the corresponding phase shift. We write 



-^[Bo + BMs) + B2w{sfj 



4(M^' + Z\2)-s 



5^/2 < 1.05 GeV, (A.5.1a) 



with A a free parameter and 



'w{s) — 



Vs - V-sq - s 



1/2 



= 1450 MeV . 



Moreover, we impose the value for the scattering length that follows from the Froissart-Gribov repre- 
sentation. We find 



Bo = (2.4±0.5) X 10^ Bi = (7.8±1.0) x 10^ 
For the inelasticity, above 1.05 GeV, 

this is negligible up to 1.25 GeV. 



B2 = (23.7±4.2) X 10^ 196±25 MeV . (A.5.1b) 

s^/^ = 1.05 GeV, e = 0.2 ± 0.2; (A.5.2c) 



A.6. The F wave 

For the F wave below s-^/^ — 1.42 GeV we fit the experimental phase shifts plus the scattering length 
as given by the Froissart-Gribov representation. We have 



cot S3 (s) — 



,1/2 



+ B0 + B: 



,1/2 



= 1.45 GeV; 



(A.6.1) 



Ba = (1.09 ± 0.03) x 10^ Si = (1.41 ± 0.04) x 10^ A = 0.051 x 10^ 



We neglect the inelasticity of the F wave below 1.45 GeV. The contribution of the F wave to all our sum 

rules is very small (but not always negligible); the interest of calculating it lies in that it provides a test 
(by its very smallness) of the convergence of the partial wave expansions. 
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A.7. The G waves 

For the GO wave, we take its imaginary part to be given by 

Im is) = BR- • 



k{Ml) J {s- Mlf + Mlr-^{k{s)/k{Ml)Y^ ' (A.7.1) 
BR = 0.17 ±0.02, = 2025 ± 8 MeV, P = 194 ± 13 MeV; c = 9.23 ± 0.46. 

For the wave G2, we can write, neglecting its eventual inelasticity, 

cot 5^l\s) = ""^'^^ B, B = (-9.1±3.3) X 10^; s^/^ > 1 GeV . (A.7.2) 

It should be noted that the expressions for the GO, G2 waves, are little more than order of 
magnitude estimates. Moreover, at low energies the expression for G2 certainly fails; below 1 GeV, an 
expression in terms of the scattering length approximation, with 

af ^ = (4.5 ± 0.2) X 10-^ M'^, 

is more appropriate. 

Appendix B: fits up to 1.42 GeV, improved with dispersion relations (CFD Set) 



In this Appendix we collect the best values for the parametrizations of the various partial waves, after 
improving with the help of dispersion relations: forward dispersion relations and Roy equations. In 
addition, we required verification (within errors) of the two crossing sum rules in Sect. 4. 

All the formulas are as before improvement, i.e., as in Appendix A; only the central values of the 
parameters change. We give in some detail only the SO and S2 waves, because now we allow variation 
of the location of the Adler zeros. For the other waves, all the formulas are exactly as in Appendix A. 

B.l. The SO wave 



The region s^/^ < 932 GeV. We take sq — 4m|-, and impose the Adler zero at s = ^Zq, with zq free. 
We find 

2\ 



cot (5q°' {s 



2k s 



iz2 



+ Bo + B- 



Vs - \/so - s 
Vs + V^o - s 



Bo 



.-\/s + V^o - s_ 



Bo =4.41 ±0.30, Bi = -26.25 ±0.60, = -15.8 ± 1.4; = 166.1 ± 4.2 MeV. 

(B.l. la) 

This corresponds to neglect of isospin breaking. If isospin breaking is taken into account as in 
ref. 20, we find instead 

Bo = 3.93 ± 0.34, Bi = -26.84 ±0.78, Ba = -10.5 ± 1.8; zo = 150.2 ± 4.2 MeV. (B.l. lb) 

The SO wave between 932 MeV and 1420 MeV. The K-matrix parameters are almost unchanged; we 
have now 



ai =0.843 ±0.017, = 0.20 ± 0.06, /3i = 1.02 ± 0.02; /Ja = 1.33 ± 0.013, 
711 =3.10 ±0.11, 712 = 1.81 ±0.05, 722 = -7.00 ± 0.04; 
Ml = 0.888 ± 0.004 GeV, M2 = 1.327 ± 0.004 GeV; 



(B.1.2) 
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B.2. The S2 wave 

The region s^/^ < 992 MeV. We here leave Z2 free. Then, 

- £^ / ^ o v/i-v/57^1 1/2 



Bo = - 80.2 ± 2.8, Bi = -69.4 ± 10.5; Z2 = 145.0 ± 3.6 MeV. (B.2.1b) 
The S2 wave between 932 MeV and 1420 MeV. We write 



cot S^^\s) = ^ _ {-^M + -Bfti [Z«ft(s) - Wh{sM)] + Bh2 [Wh{s) - Wh{sM)f^ , 

> 932 MeV; (B.2.2a) 

dwi{s) 



Bho —Bq + Biwi{sm), Bhi — Bi 



dwh{s) 
Z2 is given in (B.2.1a); 

wi{s) = ^"^{fjZI ; SI = (1050 MeV)^ 
wu{s) = ^"^^fiH ; 5, = (1450 MeV)2. 

We get 

Bfc2 = 120 ± 38. (B.2.2b) 

The inelasticity is now 

r]f\s)^l-e{l-s/sf''^, e = 0.18±0.12 (5^2 = 1.05 GeV). (B.2.2c) 

B.3. The P wave 

The region s^/^ < 2mK GeV. We have now 

Bo =1-052 ±0.011, Bi =0.17 ±0.05, Mp = 773.6 ± 0.9 MeV. (B.3.1b) 
The P wave between 2mK and 1420 MeV. The parameters are now 



Ao = 2.684 ± 0.009, Ai = 1.50 ± 0.18, A2 = -1.97 ± 0.49; 
ei =0.09 ±0.06, 62 = 0.12 ±0.11. 

B.4. Parametrization of the DO wave 

The parameters of this wave do not differ appreciably from those in Appendix A. 



(B.3.2b) 
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B.5. Parametrization of the D2 wave 

We have now 

Bo = (3.1±0.5) X 10^ 5i = (7.9±1.0) X 10^ B2 = (24.7 ±4.2) x 10^ zi = 205±25 MeV (B.5.1) 
and 

6 = 0.15 ±0.2. (B.5.2) 

B.6. The F wave 

The parameters of this wave do not differ appreciably from those in Appendix A. 
B.7. The G waves 

We have not varied the parameters of the G waves, which therefore are as in Appendix A above. 
B.8. Regge parameters 

We here give the Regge parameters, obtained with the constrained fits. They are to be used with the 
formulas of Sect. 3. 

Isospin 0. 

= 2.54 ± 0.04; cp = 0.53 ± 1.0 GeV"^; = 0.20 ± 0.10 GeV-^ (B.S.la) 

/?P' = 0.83 ± 0.05; cp' = -0.38 ± 0.4 GeV"^; ap- (0) = 0.54 ± 0.02; ap, = 0.90 GeV"^ (B. 8.1b) 
Isospin 1. 

Pp =1.30 ± 0.14; ap{0) = 0.46 ± 0.02; a'p = 0.90 GeV"^; a'^ = -0.3 GeV"'^; 

dp =2.4 ±0.5; = 0.0 ± 2.5 GeV^''. (B.8.2) 

Isospin 2. 

P2 = 0.22 ± 0.2. (B.8.3) 
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Appendix C: The conformal mapping method 

In this Appendix we explain a few of the features of the method of conformal mapping expansion. 
Although it is a standard method for calculation analytic functions with cuts (for example, it is one 
of the methods used by computers to evaluate logarithms), we hope that devoting a few lines to the 
matter would not be a waste. To make the discussion more adapted to our case, we will exemplify our 
discussion with the SO wave. 

The key point in the method is the remark that the analyticity and unitarity properties of a tttt 
partial wave amplitude/^ /(s), imply analyticity of the effective range function, tl}{s), given by (for the 
SO wave) 

=^J^ + <°-^> 

in the full complex s plane cut from —oo to 0, and from sq — Amj^ to +cxd; we are neglecting here 
inelasticity below the KK threshold. The function ip{s) is so constructed that it does not have the 
elastic cut. To find an expansion that respects this analyticity of t/j. we map the cut plane into a circle 
(Fig. C.l), which is accomplished in our case by the change of variable [conformal mapping) 




(C.2) 



w 

s 



Figure C.l. The mapping s w. 



'To lighten the notation we will, in the present Appendix, suppress indices. The rest of the notation is as in 
the main text. 
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FlGURE C.2. The w disk, \w\ < 1. The dashed Hne is the line \w\ = 0.56. 
The thick lines are the regions where one has reliable experimental data (for 
the SO wave). Images of the left hand cut and of the inelastic cut are depicted. 
The location of the sigma pole is also shown. 

Under the mapping, the left hand cut is mapped into the left half unit circle, and the inelastic 

cut into the right half of the circle; see Fig. C.2. The analyticity region (cut plane) is mapped into the 
interior of the unit circle. The function ■i/;(zf ) is then analytic, in the variable w, in the unit disk: hence, 
the analyticity properties of tp are strictly equivalent to the convergence of the Taylor expansion, 

ij{w) = Bo + Biw + B2W^ + • • • , (C.3) 

in the unit disk, < 1. 

Reverting to the variable s, the expansion (C.3) becomes, for the cotangent of the phase shift, 
the expansion 

In our work we only use the expansion in a small region, < 0.56, which is away from the 
cuts (this will be discussed in more detail later); but it is not difficult to prove that (C.4) aiso represents 
the function cot S{s) on the cuts. We show this next. 

We write a dispersion relation for the function ip^^-{s), taken to be the exact function: 

n J_oo s' - s J^^ s' -s 

^^Other methods use mapping into an ellipse, and expansions in orthogonal polynomials; see, e.g., refs. 21. 
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we neglect eventual subtractions, that play no role here. It is convenient to rewrite this with a change 
of integration variable, y — so/(2s — sq) so that we have 

, 2go f+\^ 1 Imr'^-js'jy)) 

W (s) = / ^ — . (C.5) 

TT 7_i [2s'{y) - soV s'{y) - s 

In this new variable, the left hand cut is transformed in the interval [—1,0], and the inelastic cut in 
[0,+!]. 

Next, it is clear from the presence of the functions y/s, y/s — sq that cot(5(s), or t/j{s), have the 
correct left and right hand cuts. We write a dispersion relation for the function expanded to order A'', 
ip^{s) — BnUj"", since it has cuts in the same places as the exact function, it will read 

^-W^ji^fu. ,^,, ' ■"f'-'W . (0.6, 

71" J-1 [2s'(?/) - soj^ s'{y) - s 

Now, we identify ip^ (s) and 1JJ'^^^{s) at a set of u experimental points, sj (we neglect experimental errors), 
comprised in the interval [4M^, sq]: '/'^(•Sj) ^''^^■(•Sj), j = I-, - ■ - v- This means that the functions on 
the cut, have the same averages with the set of functions 

''^^'^=[2s'iy)-so?s'iy)-s/- ^^"'^ 
y"'^'d2/^,(2/){lm^^(s'(2/))-Im^--(s'(y))}=0, j = l,...,z/. (C.8) 

Hence, the functions built with the tjj^ converge to the tp'^^', corresponding to the exact partial wave 
amplitude, both on the left hand cut and on the right hand (inelastic) cut, in the mean. In the limit in 
which one had an infinite number of experimental points, the function would be represented exactly, 
because the set of functions in (C.7) form a complete set in the interval [—1, +1] (for a proof of this in a 
physical context, cf. ref. 22). In our case we have only a finite number of points (31 for the SO case) and, 
moreover, they have experimental errors, so the representation on the cuts is valid only in the mean, 
and up to experimental errors. 

In our applications, however, this convergence on the cuts is irrelevant, as we are only fitting 
experimental data, which are located in a region away both from the left hand cut and the inelastic cut: 
cf. Fig. C.2, where we represent the experimental data we are fitting in the case of the SO wave (to the 
left we have the data obtained from ife4 decays, and if27r decays, and, to the right, those higher energy 
TTTT scattering data points included in our fits). All of them fall inside the circle \w\ < 0.56. 
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